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Abstract 

Quantum transport properties through single polycyclic hydrocarbon molecules attached to two metallic 
electrodes are studied by the use of Green's function technique. A parametric approach based on the 
tight-binding model is introduced to investigate the electronic transport through such molecular bridge 
systems. The transport properties are discussed in the aspects of (a) molecule-to-electrode coupling- 
strength and (b) quantum interference effect. 
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1 Introduction 

Molecular electronics is an essential technological 
concept of fast-growing interest since molecules con- 
stitute promising building blocks for future genera- 
tion of electronic devices where the electron trans- 
port is predominantly coherent [TJ [5] . Understand- 
ing of the fundamental processes of electron con- 
duction through individual molecules is a most im- 
portant requirement for the purposeful design of 
molecules for electronic functionalities. Electronic 
transport through molecules was first studied theo- 
retically in 1974 [3]. Following this pioneering work, 
numerous experiments [4] [5] [6] [7] H] have been per- 
formed through molecules placed between two elec- 
trodes with few nanometer separation. It is very es- 
sential to control electron conduction through such 
molecular electronic devices and the present under- 
standing about it is quite limited. For example, 
it is not very clear how the molecular transport 
is affected by the structure of the molecule itself 
or by the nature of its coupling to the electrodes. 
To design molecular electronic devices with specific 
properties, structure-conductance relationships are 
needed and in a very recent work Ernzcrhof et al. [9] 
have presented a general design principle and per- 
formed several model calculations to demonstrate 
the concept. The operation of such two-terminal 
devices is due to an applied bias. Current pass- 
ing across the junction is strongly nonlinear func- 
tion of applied bias voltage and its detailed de- 
scription is a very complex problem. The com- 
plete knowledge of the conduction mechanism in 
this scale is not well understood even today. Elec- 
tronic transport through these systems are associ- 
ated with some quantum effects, like as quantiza- 
tion of energy levels and quantum interference of 
electron waves [TQl HU H21 [T3l H21 [15l HH [17] . Here 
we focus on single molecular transport that are cur- 
rently the subject of substantial experimental, the- 
oretical and technological interest. These molecular 
systems can act as gates, switches, or transport ele- 
ments, providing new molecular functions that need 
to be well characterized and understood. 

In the present article, we reproduce an analytic 
approach based on the tight-binding framework to 
investigate the electron transport properties for the 
model of single polycyclic hydrocarbon molecules. 
Though several ab initio methods are used for the 
calculation of conductance [T51 \W\ I2UI |2"TII221 |2"31 IM1 
[25l [26] , yet simple parametric approaches [27j [28] 
129] 130] l3T| [32] are still needed for this calculation. 
The parametric study is much more flexible than 



that of the ab initio theories since the later theories 
are computationally very expensive and here we do 
attention only on the qualitative effects rather than 
the quantitative ones. This is why we restrict our 
calculations on the simple analytical formulation of 
the transport problem. 

We organize the paper as follows. With a brief 
introduction (Section 1), in Section 2, we provide a 
theoretical description for the calculation of trans- 
mission probability and current through a finite size 
conducting system sandwiched between two metal- 
lic electrodes by the use of Green's function tech- 
nique. In Section 3, we investigate the behavior of 
conductance-energy and current-voltage character- 
istics for the model of single polycyclic hydrocar- 
bon molecules and study the effects of molecule-to- 
electrode coupling strength and quantum interfer- 
ence on the above mentioned characteristics. These 
two factors i.e., the quantum interference and the 
coupling strength play crucial role on quantum 
transport through single molecular systems. Fi- 
nally, we summarize our results in Section 4. 

2 A glimpse onto the theoret- 
ical formulation 

This section follows a brief theoretical description 
for the calculation of transmission probability (T), 
conductance (g) and current (J) through a finite 
size conducting system attached to two semi-infinite 
one-dimensional (ID) metallic electrodes by the use 
of Green's function formalism. 

Let us first consider a one-dimensional conduc- 
tor with N number of atomic sites (array of filled 
black circles) connected to two semi-infinite metal- 
lic electrodes, namely, source and drain, as shown 
in Fig. [1] The conducting system within the two 
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Figure 1: Schematic view of a one-dimensional con- 
ductor with N number of sites (filled black cir- 
cles) attached to two electrodes, namely, source and 
drain. The first and the last sites of the conductor 
are labeled by 1 and N, respectively. 

electrodes can be an array of few quantum dots, or 
a single molecule, or an array of few molecules, etc. 
At much low temperatures and bias voltages, the 
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conductance g of the conductor can be written by 
using the Landaucr conductance formula. 
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where, T is the transmission probability of an elec- 
tron through the conductor. It (T) can be expressed 
in terms of the Green's function of the conduc- 
tor and its coupling to the side attached electrodes 
through the expression, 



T = Tr [T s G r c T D G a c ] 



(2) 



where, G r c and G" are the retarded and advanced 
Green's functions of the conductor, respectively. T$ 
and Tjj are the coupling matrices due to the cou- 
pling of the conductor to the source and drain, re- 
spectively. For the complete system i.e., the con- 
ductor and two electrodes, the Green's function is 
defined as, 

G=(E-H)~ 1 (3) 

where, E is the injecting energy of the source elec- 
tron. Evaluation of this Green's function requires 
the inversion of an infinite matrix as the system 
consists of the finite size conductor and the two 
semi-infinite ID electrodes, which is really a very 
difficult task. However, the entire system can be 
partitioned into sub-matrices corresponding to the 
individual sub-systems, and the Green's function 
for the conductor can be effectively written as, 



G c — (E — H c — £5 — Su) 



(4) 



where, H c is the Hamiltonian of the conductor. 
Withing the non-interacting picture, the tight- 
binding Hamiltonian of the conductor can be ex- 
pressed as, 
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where, cj (cj) is the creation (annihilation) operator 
of an electron at site i, e^'s are the site energies 
and t is the nearest-neighbor hopping integral. In 
Eq. |H E s = h\ c g s h Sc and H D = h Dc g D h' t Dc arc 
the self-energy operators due to the two electrodes, 
where gs and go are the Green's functions for the 
source and drain, respectively. hg c and h£> c are the 
coupling matrices and they will be non-zero only 
for the adjacent points in the conductor, 1 and N 
as shown in Fig. [TJ and the electrodes respectively. 
The coupling terms Y$ and for the conductor 
can be calculated through the expression, 
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where, s D ^ and £"<y D y are the retarded and ad- 
vanced self-energies, respectively, and they are con- 
jugate to each other. Datta et al. [33J have shown 
that the self-energies can be expressed like, 



^\s,D} = A {s,D} - iA{ S ,D} 



(7) 



where, A/gm are the real parts of the self-energies 
which correspond to the shift of the energy eigen- 
values of the conductor and the imaginary parts 
A{s,d} °f the self-energies represent the broadening 
of the energy levels. Since this broadening is much 
larger than the thermal broadening we restrict our 
all calculations only at absolute zero temperature. 
The real and imaginary parts of the self-energies 
can be determined in terms of the hopping inte- 
gral (i~{s,D}) between the boundary sites (1 and N) 
of the conductor and electrodes, energy (E) of the 
transmitting electron and hopping strength (v) be- 
tween nearest-neighbor sites of the electrodes. 

The coupling terms T$ and r^j can be written in 
terms of the retarded self-energy as, 
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Here, all the information regarding the conduc- 
tor to electrode coupling are included into the 
two self energies as stated above and are analyzed 
through the use of Newns- Anderson chemisorption 
theory [571 12H] ■ The detailed description of this the- 
ory is available in these two references. 

Thus, by calculating the self-energies, the cou- 
pling terms Tg and Td can be easily obtained and 
then the transmission probability (T) will be calcu- 
lated from the expression as mentioned in Eq. [5] 

Since the coupling matrices hs c and hu c are non- 
zero only for the adjacent points in the conductor, 
1 and N as shown in Fig. [TJ the transmission prob- 
ability becomes, 



T(E)=4A s n (E) A° N (E) \G 1N (E)\ 2 



(9) 



where, An =< 1|A|1 >, A NN =< N\A\N > and 
Gin =< l|G c |iV>. 

The current passing through the conductor is de- 
picted as a single-electron scattering process be- 
tween the two reservoirs of charge carriers. The 
current-voltage relation is evaluated from the fol- 
lowing expression [34], 



E F +eV/2 



I(Y) = ± J T(E) dE 



(10) 
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where, Ep is the equilibrium Fermi energy. For the 
sake of simplicity, here we assume that the entire 
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voltage is dropped across the conductor-electrode 
interfaces and it doesn't greatly affect the qual- 
itative aspects of the I-V characteristics. This 
assumption is based on the fact that the electric 
field inside the conductor, especially for short con- 
ductors, seems to have a minimal effect on the 
conductance-voltage characteristics. On the other 
hand, for quite larger conductors and higher bias 
voltages, the electric field inside the conductor may 
play a more significant role depending on the in- 
ternal structure of the conductor [33], though the 
effect becomes too small. Using the expression of 
T(E) (Eq. [9]), the final form of I(V) will be, 

E F +eV/2 

I(V) = ^ J A s n (E) A° N (E) 



E F -eV/2 
X \G 1N (E)\ 2 dE 



(11) 



Eqs.[T|l9land[TT1are the final working formule for the 
calculation of conductance g, transmission probabil- 
ity T, and current 7, respectively, through any fi- 
nite size conductor sandwiched between two metal- 
lic electrodes. 

In this article, we will study the behavior of con- 
ductance and current for some specific models of 
single polycyclic hydrocarbon molecules. Through- 
out our calculation, we use the units c = h = e = 1 . 

3 Results and their interpre- 
tation 

Here we describe conductance-energy and current- 
voltage characteristics of different single polycyclic 
hydrocarbon molecules and study the dependence 
of these characteristics on molecule-to-electrode 
coupling strength and quantum interference ef- 
fects. The schematic representations of the dif- 
ferent molecules taken into account are shown in 
Fig. [2j These molecules are: benzene (one ring), 
napthalene (two rings), anthracene (three rings) 
and tetracene (four rings). They are attached to 
the electrodes by thiol (S-H) groups. In actual 
experiments, gold (Au) electrodes are used and 
the molecules attached with them via thiol (S-H) 
groups in the chemisorption technique where hydro- 
gen (H) atoms remove and sulfur (S) atoms reside. 
In order to reveal the quantum interference effects 
on electron transport, here we consider two differ- 
ent types of bridge configurations. For one type, 
the molecules are attached to the electrodes at a-a 
sites (see the first column of Fig. [2]), the so-called 



cis configuration. In the other type, the molecules 
arc coupled to the electrodes at (3-(3 sites (see the 
second column of Fig. [5]), the so-called trans config- 
uration. 

Here we shall describe all the essential features 
of electron transport for the two distinct regimes. 
One is T/gm << t, called the weak-coupling limit 




Figure 2: Schematic representation of four differ- 
ent polycyclic hydrocarbon molecules: benzene (one 
ring), napthalene (two rings), anthracene (three 
rings) and tetracene (four rings). The molecules are 
connected to the two electrodes, at a-a positions 
the so-called cis configuration and f3-(3 positions the 
so-called trans configuration, via thiol (S-H) groups. 

and the other one is Tsg t m ~ t, called the strong- 
coupling limit, ts and td are the hopping strengths 
of the molecule to the source and drain, respec- 
tively. Throughout this work, the common set of 
values of the parameters used in our calculations 
for these two limiting cases are: t$ = To = 0.5, 
t = 2.5 (weak-coupling) and t$ = Td = 2, t = 2.5 
(strong-coupling) . The hopping integral in the elec- 
trodes is taken as v = 4 and the equilibrium Fermi 
energy Ep is fixed at 0. 

In Fig. [31 we plot the conductance (g) func- 
tion of the injecting electron energy (E) for the 
molecular bridge systems where the molecules are 
coupled to the electrodes in the trans configura- 
tion, (a), (b), (c) and (d) correspond to the re- 
sults for the benzene, napthalene, anthracene and 
tetracene molecules, respectively The solid curves 
denote the results for the weak molecule-electrode 
coupling case and it is observed that the conduc- 
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tance shows sharp resonant peaks for some par- 
ticular energy values, while it (g) drops to zero 
almost for all other energies. At resonance, the 
conductance approaches to 2 so that the transmis- 
sion probability (T) becomes unity (from the Lan- 
dauer conductance formula g = 2T, see Eq. [T]with 
e = h = 1). The resonant peaks in the conductance 
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conductance-energy characteristics for the molec- 
ular bridge systems where the molecules are con- 
nected to the electrodes in the cis configuration, 
(a), (b), (c) and (d) correspond to the results for 
the benzene, napthalenc, anthracene and tetracene 
molecules, respectively, where the solid and dot- 
ted curves indicate the same meaning as in Fig. [3l 
These results show that some of the conductance 
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Figure 3: g-E spectra for the polycyclic hydrocar- 
bon molecules connected to the electrodes in the 
trans configuration, (a), (b), (c) and (d) corre- 
spond to the results for the benzene, napthalenc, 
anthracene and tetracene molecules, respectively. 
The solid and dotted curves represent the results 
for the weak and strong molecule-to-electrode cou- 
pling limits, respectively. 

spectrum coincide with the eigenenergics of the sin- 
gle hydrocarbon molecules. Therefore, the conduc- 
tance spectrum manifests itself the electronic struc- 
ture of the molecules. On the other hand, in the 
strong-coupling limit the resonances get substantial 
widths, as shown by the dotted curves in Fig. [3] and 
it emphasizes that the electron conduction takes 
place almost for all energy values. This is due to 
the broadening of the molecular energy levels, where 
the contribution comes from the imaginary parts of 
the self-energies Egm) [33] . as mentioned earlier in 
the strong-coupling case. 

To characterize the quantum interference effects 
on the electron transport, in Fig. [4j we plot the 



Figure 4: g-E curves for the polycyclic hydrocar- 
bon molecules connected to the electrodes in the 
cis configuration, (a), (b), (c) and (d) correspond to 
the results for the benzene, napthalenc, anthracene 
and tetracene molecules, respectively. The solid 
and dotted lines represent the similar meaning as 
in Fig. H 

peaks do not reach to unity anymore and get much 
reduced value. This behavior can be understood as 
follows. The electrons are carried from the source to 
drain through the molecules and thus the electron 
waves propagating along the two arms of the molec- 
ular ring(s) may suffer a phase shift between them- 
selves, according to the result of quantum interfer- 
ence between the various pathways that the electron 
can take. Therefore, the probability amplitude of 
the electron across the molecules becomes strength- 
ened or weakened. It emphasizes itself especially as 
transmittance cancellations and anti-resonances in 
the transmission (conductance) spectrum. Thus, it 
can be emphasized that the electron transmission 
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is strongly affected by the quantum interference ef- 
fects, and hence, the molecule-to-electrode interface 
structures. 

The scenario of electron transfer through the 
molecular junction is much more clearly observed 
from the current-voltage characteristics. Current 
through the molecular systems is computed by the 
integration procedure of the transmission function 



staircase-like structure with sharp steps. This is 
due to the discreteness of molecular resonances as 
shown by the solid curves in Fig. [3J As the volt- 
age increases, the electrochemical potentials on the 
electrodes are shifted and eventually cross one of 
the molecular energy levels. Accordingly, a current 
channel is opened up and a jump in the I-V curve 
appears. The shape and height of these current 
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Figure 5: I-V characteristics for the polycyclic hy- 
drocarbon molecules connected to the electrodes in 
the trans configuration, (a) and (b) correspond to 
the results for the weak and strong molecule-to- 
electrode coupling cases, respectively. The solid, 
dotted, dashed and dot-dashed lines represent the 
currents for the benzene, napthalcnc, anthracene 
and tctraccnc molecules, respectively. 

T (see Eq. [TTj) . The behavior of the transmission 
function is similar to that of the conductance spec- 
trum since the equation g = 2T is satisfied from 
the Landauer conductance formula. In Fig. we 
plot the I-V characteristics for the hydrocarbon 
molecules connected to the electrodes in the trans 
configuration, (a) and (b) correspond to the cur- 
rents for the molecular bridge systems in the weak- 
and strong-coupling limits, respectively. The solid, 
dotted, dashed and dot-dashed lines represent the 
results for the benzene, napthalcnc, anthracene and 
tetracene molecules, respectively. It is observed 
that, in the weak-coupling case the current shows 



Figure 6: I-V spectra for the polycyclic hydrocar- 
bon molecules connected to the electrodes in the cis 
configuration, (a) and (b) correspond to the results 
for the weak and strong molecule-to-electrode cou- 
pling cases, respectively. The solid, dotted, dashed 
and dot-dashed lines represent the identical mean- 
ing as in Fig. [5] 

steps depend on the width of the molecular reso- 
nances. With the increase of molecule-to-electrode 
coupling strength, the current varies almost contin- 
uously with the applied bias voltage and achieves 
much higher values, as shown in Fig. E^b). This 
continuous variation of the current is due to the 
broadening of the conductance resonant peaks (see 
the dotted curves of Fig. [3]) in the strong molecule- 
to-electrode coupling limit . 

The quantum interference effects on the electron 
transmission through the molecular bridge systems 
is much more clearly visible from Fig. [6l where 
we plot the I-V characteristics of the hydrocarbon 
molecules connected to the electrodes in the cis 
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configuration, (a) and (b) correspond to the cur- 
rents in the two limiting cases, respectively. The 
solid, dotted, dashed and dot-dashed curves give 
the same meaning as in Fig. [5j These results predict 
that the current amplitudes get reduced enormously 
compared to the results predicted for the molecules 
connected to the electrodes in the trans configu- 
ration. This is solely due to the quantum interfer- 
ence effects between the different pathways that the 
electron can take. Therefore, we can predict that 
designing a molecular device is significantly influ- 
enced by the quantum interference effects i.e., the 
molcculc-to-electrodc interface structures. 



4 Concluding remarks 

To summarize, we have introduced a parametric ap- 
proach based on the tight-binding model to investi- 
gate the electron transport properties of four differ- 
ent polycyclic hydrocarbon molecules sandwiched 
between two metallic electrodes. This technique 
can be used to study the electronic transport in any 
complicated molecular bridge system. Electronic 
conduction through the hydrocarbon molecules is 
strongly influenced by the molecule-to-electrode 
coupling strength and the quantum interference ef- 
fects. Our study provides that designing a whole 
system that includes not only the molecule but also 
the molcculc-to-clcctrodc coupling and the inter- 
face structures are highly important in fabricating 
molecular electronic devices. 

More studies are expected to take the Schottky 
effect, comes from the charge transfer across the 
metal-molecule interfaces, the static Stark effect, 
which is taken into account for the modification of 
the electronic structure of the molecular bridge due 
to the applied bias voltage (essential especially for 
higher voltages). However, all these effects can be 
included into our framework by a simple general- 
ization of the presented formalism. Here, we have 
also neglected the effects of inelastic scattering pro- 
cesses and electron-electron correlation to charac- 
terize the electronic transport through such molec- 
ular bridges. 
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